Draft version May 9, 2008 

Preprint typeset using I£T^X style cmulateapj v. 10/09/06 



OPTIMAL COADDITION OF IMAGING DATA FOR RAPIDLY FADING GAMMA-RAY BURST 

AFTERGLOWS 

A. N. Morgan 1,2 , D. E. Vanden Berk 1 , P. W. A. Roming 1 , J. A. Nousek 1 , T. S. Koch 1 , A. A. Breeveld 3 , M. de 
Pasquale 3 , S. T. Holland 4 ' 5 , N. P. M. Kuin 3 , M. J. Page 3 , M. Still 3 

Draft version May 9, 2008 

ABSTRACT 

We present a technique for optimal coaddition of image data for rapidly varying sources, with specific 
application to gamma-ray burst (GRB) afterglows. Unweighted coaddition of rapidly fading afterglow 
lightcurve data becomes counterproductive relatively quickly. It is better to stop coaddition of the 
data once noise dominates late exposures. A better alternative is to optimally weight each exposure to 
maximize the signal-to-noise ratio (S/N) of the final coadded image data. By using information about 
GRB lightcurves and image noise characteristics, optimal image coaddition increases the probability 
of afterglow detection and places the most stringent upper limits on non-detections. For a temporal 
power law flux decay typical of GRB afterglows, optimal coaddition has the greatest potential to 
improve the S/N of afterglow imaging data (relative to unweighted coaddition), when the decay 
rate is high, the source count rate is low, and the background rate is high. The optimal coaddition 
technique is demonstrated with applications to Swift Ultraviolet/Optical Telescope (UVOT) data of 
several GRBs, with and without detected afterglows. 
Subject headings: methods: data analysis; gamma-rays: bursts 



1. INTRODUCTION 

The early detection of gamma-ray burst (GRB) op tical afterglows is c rucial to fulfill several key science goals set 
forth by the Swift Ga mma-ray Burst Explorer mission (|Gehrels et al.ll2004h . The Ultraviolet /Optical telescope (UVOT; 
iRoming et alj [2005a! ) onboard Swift was designed to capture these early (~ 1 minute after burst trigger) afterglows 
in order to both study the behavior of the early optical afterglow and to provide accurate sub-arcsecond positions 
for ground-based follow-up observations. How ever, approximately 60% of the GRBs detected by Swift lack an optical 
counterpart in rapidly available UVOT data (|Marshal]||2007l : IRoming fe Masonl 120061 : iMason et ai1l2005D . and more 
than 40% lack an optical or IR counterpart in data from any telescope. This unexpectedly high fraction of optical 
non-detections highlights the importance of ensuring that the available data are being used to their full potential. 

In particular, it is imperative that data are coadded in such a way that the chance of revealing a faint detection 
is maximized or the most stringent upper-limits are obtained. The rapidly fading behavior of GRBs ensures that 
coadding more data will eventually become detrimental as the burst signal fades below the backgro und level. Upper 
limit s have often been reported via the Gamma-ray bursts Coordinates Network (GCN) Circulars (jBarthelmv et al.l 
fl9Mli"99i which, in our experience, are often calcula ted using the coaddition of an arbitrary number of unweighted 
images (e.g.. iMorgan et aT1 l2005: Ro ming et al.l [200 5b). By continuing to coadd capriciously, we run the risk of not 
reporting the most useful information to the rest of the GRB community and may even overlook a faint afterglow. In 
order to best use the imaging data from Swift, a method for determining when to stop coadding burst data, or for 
optimally coadding the data, must be implemented. 

This paper describes the necessary criteria and the technique for optimally coadding faint GRB afterglow imaging 
data, and for improving detection limits when no afterglow has been detected in the individual images. Equations 
are derived that describe the signal-to-noise ratio (S/N) of a number of summed exposures for an object whose flux 
density fades according to a simple power law decay 

F y oc t a , (1) 

where a < is the temporal decay index. The S/N of an optimally coadded imaging data set depends on the decay 
index, the start and stop times of each exposure, the measured counts in each exposure, and the noise levels of each 
exposure. For bursts which lack an optical afterglow detection, however, the value of the crucial decay index parameter 

is unknown, making it necessary to adopt assumptio ns in this case. 

Observations with SwijVs X-ray Telescope (XRT; iBurrows et al.l [20051) suggest a canonical GRB X-ray afterglow 
lightcurve, consisting of three distinct power law segments with decay indices ranging from very steep decay at early 
times (-5 < a < -3;t < 300 s), followed by a shallow decay (-0.8 < a < -0.2) unt il i ~ 10 4 — 10 5 s, and finally 
a decay of medium (—1.5 < a < —1) steepness ([Nousek et al.l [20061 : IZhang et al.l l2006f) . Optical and UV afterglows 
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typically display an initial phase lasting up to about 500 s during which the lightcurve is slowly deca ying or even 
rising , and a longer steeper power law decay phase with a temporal slope of about a = —0.9 ± 0.4 (jOates et al.l 
2008) . Flaring at early tim es, as is often seen in X-ray lightcurves (e.g . Nousek et al. | l2006ft , is occasionally manifested 



in op tical lightcurves (e.g. [Holland et ail 120021 ; Uakobsson et al.l 12004 ; iBlustin et aTl l2006; W ei et al.|[200l iDai et all 



2007). Optical and X-ray lightcurves may, but typically do not, track each other closely. We assume the simplest 
temporal decay model given by equation ([1]) in our modeling of S/N, but it is straightforward to generalize the 
technique to arbitrary lightcurves. 

In § ??, observations of GRB afterglows using the Swift UVOT are described. The optimal image data coaddition 
technique is developed in § ??. The technique is tested with simulated data in § ??, and applied to examples of real 
UVOT afterglow data in § ??. A discussion and summary are given in § ??. 

2. SWIFT UVOT OBSERVATIONS OF GRB AFTERGLOWS 

The techniques described here apply generically to observations made by a wide variety of telescopes and detectors, 
however, the UVOT telescope on board the Swift satellit e makes rapid and long -term observations of almost all GRB 
fields identified by the Swift Burst Alert Telescope (BAT; lBarthelmv et al.| [2005). For this reason the specific examples 



used here will focus on data obtained with the UVOT. In this section we briefly describe the UVOT and the GRB 
afterglow data sets it typically produces. 

The Swift UVOT is a 30 cm diameter telescope with a 17 x 17 arcminute field of view, with a detector sensitive to wave- 
lengths bet ween approximately 1 600 and 8000A, covered by seven photometric filters labeled u yw2, uvm2 , uyml . u, b, v, 
and white (Roming et al. 2005a). The photometric calibration of the UVOT is described bv lPoole et alJ (|2008h . 

The average time to the first UVOT observation of a GRB location after the BAT trigger is 110s, barring constraints 
due to the positions of a burst relative to the Sun, Earth, and Moon. The automated GRB observing sequence utilizes 
all seven filters, with the order and exposure times determined with respect to the time of the burst. The simulations 
described in § ?? are based on a typical set of exposures in the UVOT v band, up to about 10 5 s post-burst. A normal 
set of v band exposures would consist of a 400 s "finding chart" exposure shortly after slewing to the GRB position, a 
series of short (10 to 20 s) exposures, another finding chart exposure, another series of short exposures, and finally a 
set of longer exposures (100 to 900 s). The sequence described here is typical for the majority of Swift detected GRBs 
and is currently implemented for the detection of new bursts. Observations often continue to be made for several days 
or weeks, depending on the flux of the afterglow seen by the XRT. It is very rare, however, for any significant flux to 
be detected at optical or UV wavelengths at these late times; so late-time observations will not be considered for the 
simulations in § ??. 

GRB afterglows are detected in the first white or v observations approximately 40% of the time when those observa- 
tions are made less than 500 s after the burst. An analysis of the im plications for the nature of GRB afterglows based 
on the UVOT detection statistics is given by iRoming et alJ (|2006f ) . Additional detections are sometimes made after 
the coaddition of frames taken over a period of time. In the next section, we examine the effectiveness of standard 
unweighted frame coaddition of GRB afterglow data, and describe a method for optimal coaddition of imaging data 
of temporally varying sources. 

3. OPTIMAL COADDITION OF GRB IMAGE DATA 
3.1. Unweighted Coaddition 

Assuming a source with a lightcurve that follows a power law decay, we derive an equation that gives an estimate 
of the S/N in an aperture surrounding the source for an unweighted sum of n exposures. The equation for the final 
S/N depends on the initial source count rate of the afterglow, the background count rate for each exposure, the start 
and stop times of each exposure, and the temporal decay index of the burst. With this equation, one can calculate 
the exposure at which the maximum S/N occurs for a given burst if unweighted coaddition is adopted. We implicitly 
assume that the measurements are made in the same aperture surrounding the same known source location in each 
of the images. Summing the images and then measuring counts in an aperture in the summed image is equivalent 
to summing the measurements from apertures in the individual images; here we proceed as if individual aperture 
measurements are being summed. 

There are several sources of noise for a given detector. Among them are noise due to the signal itself (N src ), noise 
due to the sky background (N s k y ), dark noise (Nd ar k), and read noise (V reQt j). Assuming all sources of noise are known 
to perfect accuracy and are uncorrelated, the total noise is the sum in quadrature of these four quantities: 

N 2 = N 2 src + N 2 sky + N 2 ark + N 2 ead . (2) 

For the UVOT, th e read noise is zero since it is a photon counting device, and the dark noise is insignificant 
(M ason et aT1l2004l ). and thus the dominant contributors to the noise are the source and the sky background. We 
assume that the source and sky counts are Poisson distributed, but that the summed image counts are large enough 
that Gaussian statistics apply, so that the total noise is simply the square root of the number of counts: 

N 2 = ^ (C src .i + C s kys) , (3) 
i 

where C src> i and C s k y ,i are the measured number of source and sky background counts in the i th exposure, respectively. 
The signal-to-noise ratio of the summed image is thus given by 

S/N sum = ■ (4) 
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The estimated total number of counts in an exposure from a given source is the integral of the count rate R(t) from 
the start to stop times of an exposure. If we approximate the sky background count rate (R s ky,i) &s constant during 
the i th exposure, the estimated sky counts in that exposure are 



Csky,i — Rsky,i (tstop.i ^start,i) 



(5) 



where t s t ar t,i and t st0 p,i are the start and stop times of the i exposure after the initial burst trigger. 

Using the assumption of a simple temporal power law decay for gamma-ray burst afterglows, the model source count 
rate Rmodei 

can be parameterized as 

Rmodel(t) = Rl ( ~ 



(6) 



where a is the temporal decay index and R\ is the initial source count rate, defined to be the average count rate of 
the I s * exposure given by 

Ri= (t C Z C t y 00 

\ L stop,l L start,lJ 

The parameter t\ is thus the weighted midpoint of the exposure at which time the count rate is R\. Integrating ([S]), 
the model number of source counts in the i th exposure is given by 
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For the remainder of the derivation, we shall assume for simplicity that a is never exactly -1. Inserting |(7J) into 
i = l and solving for t\ gives 
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Inserting © into §8§ gives 
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The time since the burst associated with the coadded data is not well defined, since it involves many different time 
intervals, each sampling a different portion of a light curve. However, a reasonable definition for a characteristic post 
burst time, t, is the count weighted average time of the observations. In other words, i is the average arrival time 
of the counts. Because the count rate is expected to change across an exposure, the observing times are weighted by 
the integrated model count rate across each exposure window, rather than weighting them with the measured counts. 
Thus we arrive at the characteristic time 
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which is independent of the initial count rate. 

If the temporal decay model is a good approximation to the real lightcurve, the S/N for a summed set of exposures 
can be predicted as 
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With this equation, one can predict the signal-to-noise for a summed image of n exposures by simply knowing the start 
and stop times of each exposure (relative to the time of burst) and measuring both the initial source count rate and 
the background count rates for each exposure. In typical sets of observations, the S/N will initially increase with the 
coaddition of more frames, but as the count rate drops the noise will increase faster than the signal, until eventually 
the S/N will reach a peak value, then begin to fall. By modeling the S/N, one can determine the amount of data to 
coadd to achieve the most significant detection for unweighted data. 
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3.2. Weighted Exposures 

While equation (fT5)) can be used to predict a peak S/N in unweighted coadded data, the significance of a source 
detection in a series of images can be raised by optimally weighting the data in each exposure. Additionally, with 
optimal weighting, the S/N will not drop as the source cou nts become le ss significant. The derivation of the optimal 
weighting scheme parallels in part the arguments given by iHornel (|1986l ) for the optimal extraction of spectroscopic 
data. In that case, the spatial profile of a spectrum at a given wavelength was used to provide independent estimates, 
at each pixel, of the total source counts in the spectrum. Each estimate was weighted in such a way that the average 
value gives the minimum variance in the the total source counts. In the case of GRBs, the temporal lightcurve takes 
the place of the spatial spectrum profile, and each observation is the synonym of a pixel in the spectrum profile. 

Let Pi be the probability that a detected photon from a GRB afterglow is recorded in frame i. The sum of the 
probabilities over all the n frames is normalized to unity 

n 

X>i = L (14) 

i=l 

The probability for frame i is given by the expected GRB afterglow counts in that frame, C mo dei,i (given, for example 
by equation (jTUJ) ) , divided by the total number of expected counts, which is the sum over the expected counts in all 
of the frames 

Pi C nlo del,i , . 

i = ^n^ ■ (15) 

/ „ • ^model.i 

The probability function is the normalized model GRB afterglow lightcurve. We assume here that the lightcurve is 
known, but we will discuss the case in which the lightcurve is uncertain in § ??. 

Given the normalized lightcurve, an estimate of the total number of source counts can be made for each frame, by 
dividing the source counts by the probability function for that frame 

Ctot.i ~ C srCy i/Pi . (16) 

The average value of all of the estimates is then a linear and unbiased estimator of the total number of source counts 
(Horne1 [l986D . Including a weighting factor Wi, for each frame, the estimator of the total source counts Ctot is 

„ _ Ei w iC S rc,i/ Pi 

'-■tot — • I 1 ' J 

The total count estimate is optimal when the variance on the estimator is minimized. The variance of equation (|17|) . 
assuming equation ([3]), is 

1 - 

V(Ctot) = — -a £ wi (C srv + C skyii ) /Pi , (18) 
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modulo a multiplicative constant, so that equation (fTT)) becomes 
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For a constant source, equation (|20[) reduces to the familiar inverse variance weighting. The optimal weighting factor 
(fT^|) not only accounts for differences in the variances of each frame, but weights each frame differently depending on 
the number of source counts expected in the frames. The S/N of the coadded weighted image data is 

q/M _ E7PiCsrc,i/(C arc , i +C sky , i ) 

*J / ly sum, weight — n 2 1/2 ' / 

[Ei Pi I (CsrcA + Cskyd)] 

In contrast to the unweighted S/N given in equation the weighted sum S/N runs little risk of decreasing as the 
count rate drops, if the lightcurve model is sufficiently accurate. Frames with few expected counts contribute little 
to either the final signal or noise. The coaddition process does not need to stop at a particular frame to achieve a 
maximum S/N - statistically, the maximum S/N will always be reached with the coaddition of all of the frames. 

The expected number of counts in each frame is found from the lightcurve model. Assuming again that GRB 
afterglow lightcurves can be described by a single power law, as in equation (|10|) . and using equation (|15|) for the 
probabilities Pi, the optimally weighted total number of afterglow counts in a set of exposures is 
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and the S/N is 
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This shows that Ctot can be determined independently of the initial model count rate R\ , since the equation instead 
contains the decay slope and observing times. The total count estimate does not correspond to a particular count 
rate, because the count rate changes over time. However, the initial estimated count rate can be found by setting the 
weighted total counts equal to the sum of the model counts in equation (|10|) . and solving for R\. Then the count rate 
at any moment can be estimated from the lightcurve model in equation (|6j). 

The characteristic time since the burst of the observations can be defined in a manner similar to equation (jTTJ) , 
except that in this case each exposure is further weighted by its optimal weight, Wi, 
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For the optimal weights defined by equation (|19p . the characteristic time is 



t ■ 



a+2 
stop.i 



t a 

start a 



)(* 



a+l 
stop,i 



i) /(C s 



Csky,i ) 



which is independent of the initial count rate. 
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3.3. Detection Limits 

When an afterglow is not detected in any single image, coaddition of the image data may improve the final S/N 
enough to provide a detection. Even when the final coadded image does not yield a detection, optimal coaddition 
can be used to place stronger detection limits on the dataset, than would result from unweighted coaddition. If 
the lightcurve shape can be assumed, the source counts C src ^, can be replaced by the lightcurve model C moc iei.i in 
equations Q and (|2"Tj) . and the equations set to a minimum S/N (e.g. 3) required for a detection, then solved for the 
parameters of the lightcurve model. For the case of a power law decay with a given index (equation (fT0|0. this will 
yield the maximum initial average count rate, R\, m ax that could have occurred without producing a detection in the 
coadded data. In other words, given a non-detection in the coadded data, R\. m ax is the limiting measured count rate 
of the source during the time of the first exposure. 

For unweighted data, the expression for R\ im ax can be derived by solving equation (|13p for R\ t max- The solution 
is quite simple in that case, since the equation is a quadratic. The unweighted maximum initial count rate may be 
either lower or higher (better or worse respectively) than the detection limit of the first image by itself, depending on 
how many unweighted noise-dominated exposures are coadded. 

When the data are optimally coadded, equation (|2ip applies, with C src ,i replaced with C mo dei,i- Assuming a power 
law decay model for the lightcurve, the S/N is described by equation (|23| : li in place of C srCtl , that 

equation cannot be solved analytically for Ri. m ax, but it is not difficult to solve numerically. Initial (but slightly 
over-estimated) guesses for R\^ m ax can be made by assuming that the source counts in the denominator of equation 
(|23|) are negligible compared to the background counts, which yields an equation that can be solved analytically for 
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where S/N m i n is the minimum S/N required for the detection of a source. By utilizing the optimally coadded data, 
the initial average count rate limit will always be more sensitive than the limit of the first image by itself. Statistically, 
the optimally weighted R\ t max will be more sensitive than the value calculated from the unweighted data. 



4. TESTS WITH SIMULATED GRB AFTERGLOW DATA 

The final coadded S/N of GRB afterglow data depends upon the observing times, the initial afterglow count rate, 
the temporal decay slope, the background count rate, and the count rates in each of the images. To test how the 
quality of coadded data depends on these parameters, we have constructed simulated GRB afterglow lightcurve data, 
and applied both unweighted and optimally weighted coaddition. While focusing on typical UVOT parameters, the 
qualitative results of the simulations apply generally to most GRB afterglow observations. 

For each simulation, the decay slope a, the initial afterglow count rate R\, and the background count rate C s k y 
(assumed constant for simplicity) were specified. The observing times were kept the same in each case, and are 
representative of a typical UVOT observing sequence in the v band (See §??). In each observing time interval, the 
lightcurve model was integrated over the exposure times to find the expected number of counts, which was used as 
a Poisson mean. The background count rate was also integrated over the exposure times to find the Poisson mean 
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background counts for each time interval. The simulated detected counts from both the afterglow and background 
were drawn from Poisson distributions, given their respective mean values. The sum of the detected afterglow and 
background counts is the total simulated observed counts. The mean background counts (assumed to be accurately 
measurable) were subtracted from the total counts to give the simulated measured afterglow counts. Using these 
values, the coadded S/N at each time interval was calculated according to equation (j4|) for the unweighted case, and 
equation (|23|) for the optimally weighted case. 

Optical temporal dec ay slopes for GRB afterglows are usually observed to be in the range a = —0.5 to —2.0, with an 
average of about -0.9 (pates et al.ll2008D . To see how the effectiveness of optimal coaddition is affected by different 
decay slopes, data were coadded for simulated lightcurves with a range of decay slopes, while the other parameters 
were held constant. Figure[Tj shows the cumulative S/N after coadding the data through each observing interval, using 
both unweighted, and optimally weighted coaddition, for two different values of the decay slope. Error bars show the 
68.3% (nominal la) confidence limits from the results of 1000 simulations. In each case, the unweighted S/N reaches 
a peak value, then decreases as more data are coadded. This happens because the noise dominates the observations 
to a greater extent with each successive observation. In contrast, the S/N of the optimally weighted data rises, then 
flattens to a nearly constant value at long post-burst times. The late-time data contributes very little signal or noise, 
so the S/N values remain almost unchanged. The effect is the same, regardless of the value of the decay slope, but 
the relative difference between the final weighted and unweighted S/N is greatest for the data with the steepest decay 
slope. When the lightcurves are steep, the afterglow fades below the noise level more quickly, so the unweighted S/N 
drops more rapidly, while the optimally weighted S/N flattens out sooner. 

Figure[2] shows how the coadded S/N is affected by differences in the initial count rate; the decay slope and back- 
ground count rates were fixed at typical values. The results show that the final weighted S/N is significantly better 
than the unweighted S/N in all cases, and that optimal coaddition has the greatest potential to improve the final 
S/N when the count rate is low. These statements are true in a statistical sense, but at low count rate, there is no 
guarantee that the cumulative weighted S/N at any observing time will exceed the unweighted value, as the dispersion 
of the simulations shows. Nonetheless, as shown by Fig.[2j in typical cases optimal coaddition is more likely to result 
in a significant (S/N > 3) detection, than when using unweighted coaddition. 

The background count rate clearly affects the detectability of fading afterglows. The background in UVOT images 
(although low compared to ground-based telescope images) can change greatly from filter to filter, and even among 
images using the same filter, depending on the directions of the Sun, Earth, and Moon. Figure[3] shows that optimal 
coaddition improves the S/N in the presence of a wide range of background rates (with decay slope and initial count 
rate fixed in the simulations), particularly when the background rate is high. When background count rates are very 
low, as is often the case with X-ray telescopes, there is little benefit from optimal coaddition. For most optical and 
infra-red detectors, the background is usually the dominant noise factor, so optimal coaddition is worth the minimal 
implementation effort. 

Often the decay slope is unknown, or is not known with great accuracy, so it is of interest to know how the final 
S/N depends upon the decay slope used in the optimal coaddition. Simulations were made using typical values of 
the true decay slope a = —1.2, the initial afterglow count rate R\ = 3 counts s _1 , and the background count rate 
C s ky — 2.5countss _1 . The model decay slope a mo dei (not necessarily the same as the true slope), used for the optimal 
coaddition was varied over a wide range of values. The median final S/N in the coadded data simulations were 
calculated for each value of a mo( i e i- The results are shown in Fig.|H As expected, the final S/N is maximized when 
the model decay slope closely matches the true slope. Figured] shows that, at least for typical cases, there is some 
leeway in the value of the model decay slope; differences of up to a few tenths in the slope change the final S/N only 
a small amount. However, the model slope value cannot be arbitrarily different from the true value, or the final S/N 
can actually drop well below the maximum S/N in the unweighted case (shown by a dashed line in Fig.[4|). 

The results shown in Fig. 3] suggest a way to estimate the decay slope when it is not known. Optimal data coaddition 
can be repeated for a data set, for a range of values of the model decay slope. The value that maximizes the final 
S/N is most likely to come closest to the true decay slope. For the example given here, we calculated the value of the 
decay slope that maximized the S/N in the optimally coadded data, for each simulated light curve. The mean value 
of the decay slope was a = —1.20 with a standard deviation among the simulations of a(a) — 0.19. The decay slope 
estimation technique works well even in cases for which an afterglow is detected only after optimal coaddition . For 
example, in simulations where the afterglow decay slope is a = —0.9 (the average given bv lOates et"ai] (|2008f n. the 
initial count rate is R\ — 0.9countss _1 , and the background count rate C s k y = 5.0 counts s , the optimally coadded 
S/N is just over 3, while the estimated value of the decay slope is a = —0.89 ± 0.44. Thus, this technique can be used 
to calculate reasonable estimates of the decay slope, even in low S/N circumstances. 

As mentioned earlier, X-ray and optic al lightcurves usually do not track each other very well, particularly at early 
times. There are exceptions t o this (e.g. iBlustin et al.|[2006t iGrupe et al.ll2007l ). but usually the X-ray decay slope is 
steeper than the optical slope (jOates et al.ll2008h . so at best the X-ray slope can be used as a lower limit for the optical 
slope estimate. When the optical slope is unknown, the best approach for optimal coaddition is probably to assume a 
range of reasonable decay slopes, bracketed at the steep end by the X-ray slope, and to calculate limiting magnitudes 
in each case. 
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5. APPLICATION TO SWIFT GRB IMAGE DATA 
5.1. Optically Bright Bursts 

To demonstrate the effect of optimal coaddition on real data, imaging data from the bright well-sampled optical 
afterglow of GRB 050525A were coadded using the weighted and unwe ighted methods. T he GRB was discovered by 
the Swift BAT, and the spacecraft slewed to the source immediately (Ban d et al.l 12005). The UVOT began obser- 
vations about 65s af ter the burst, and a bright afterglow was detected in images taken with all seven UVOT filters 
(jHolland et al.ll2005h . UVOT imaging data of the afterglow field were collected up to about 1.2 x 10 6 s after the burst; 
we have used data only up to a few x 10 5 s after the burst, since the much later time data are e ntirely consistent wit h 
no source flux and add nothing to this example. Full analysis of the Swift data was presented bv lBlustin et~aT1 (|2006l) . 
Here we show how the S/N of UVOT data from a well-detected afterglow can depend strongly on how the data are 
coadded. 

The UVOT v band lightcurve of GRB 050525A is shown in the top panel of Fig.[5l and was constructed as follows. 
Source plus background counts in the individual v band images were measured using a 3 arcsecond radius aperture. A 
large background annulus region, centered on the afterglow position, and excluding small regions around any detected 
stars near the annulus, was used to estimate the background counts. Source counts were measured by subtracting the 
expected average background counts inside the source extraction region from the total counts inside that region. A 
power law model was fit to the lightcurve by finding the parameters that minimized the x 2 value of the fit, using the 
linear count rate data and the associated uncertainties. The best fit power law has an index a = —1.03, and is shown 
in Fig. [5] The value for the power law index is fairly close to that found bv lBlustin et al.l (l2006f). who fo und a value of 
a = —1.14 when a simultaneous fit was made to data in multiple bands. However. iBlustin et al.l (2006) found that a 
single power law does not give a full statistical description of the optical/UV lightcurve in any band. For this example, 
we have ignored the deviations from a single power law, since they would unnecessarily complicate the analysis, and 
would provide very little improvement to the optimal coaddition of the data. 

The v band data were coadded without weighting, and using optimal weighting according to equation |22|) . The 
GRB afterglow counts expected in each image were calculated according to equation (fit))) , using the best fit value 
of the decay index in the v band. The cumulative S/N at each exposure was calculated from equation (|13[) for the 
unweighted data, and from equation (|23|) for the weighted data. The bottom panel of Fig. [5] shows the cumulative S/N 
for both the weighted and unweighted cases. The S/N values in both cases match each other very closely through the 
first w 1000s after the burst. Within that time range, the afterglow was still quite bright, and the noise in the images 
is dominated by the counts from the afterglow itself; optimal coaddition holds little advantage in that case. However, 
at later times, when the afterglow count rate had dropped to relatively low values, continued unweighted coaddition of 
the data becomes detrimental. When all of the data are coadded without weighting, the final S/N is well below that 
of the first image by itself. When the data are optimally coadded, the S/N at later times actually increases slightly, 
and levels off at a maximum value. 

While coaddition of data in the case of GRB 050525A was not necessary for a detection, this example illustrates 
the danger of coadding image data without considering the rapidly decaying nature of afterglow lightcurves. At best, 
unweighted coaddition makes no improvements, and at worst it can severely reduce the quality of coadded data. 
Optimal coaddition nearly guarantees that the highest S/N will be returned from coadded data, regardless of how 
many images are coadded. In the next two sections, we illustrate the benefits of optimal coaddition applied to data 
with low S/N in individual images. 

5.2. Optical Detections from Optimally Coadded Data 

A considerable fraction of optical GRB afterglows are not significantly detected (at > 3cr) in individual UVOT 
images taken with a specific filter, but may be significantly detected in the coadded data. An example in which 
the UVOT v band data reveal no afterglow in individual images, but optimal coaddition results in a detection is 
GRB 060604. This GRB was di scovered by the Swift BAT, and the spacecraft slewed to the source position about 100s 
after the detection (Page 2006). An afterglow candidate was found in the initi al UVOT white filter image, and later 
confirmed in unweighted coad ded UVOT b and u images (jBlustin fc Pagdl2006f ). and by ground-based R band images 
(jTanvir. Rol. fc Hew ctt 2006]). No significant afterglow was found in any of the individual UVOT v band images; 
unwe ighted coaddition of over 1100s of data yielded a possible source with a significance of 2.4 a (jBlustin fc Page! 
2006). This case is a good test of the optimal coaddition technique, since we know that there was an optical afterglow 
in several bands, but that the source was only marginally detected in another band. The optimal coaddition technique 
was applied to the v band data to determine if the significance of the afterglow could be improved. 

The v band data consist of 17 images, excluding the initial very short settling image (taken while the spacecraft is 
still slewing), starting about 220s after the burst, and ending about 10 5 s after the burst. The optimal weights used in 
equation (|17p were determined from measurements of the source counts and background counts in each image, and the 
temporal decay slope. Source plus background counts in the v band images were measured inside an aperture centered 
on the known GRB afterglow position, with an aperture of radius 3 arcseconds , which is about optimal for maximizing 
the S/N of faint sources in UVOT images (jPoole et al.ll2008t iLi et al.l l2006aT ) . A large background extraction region 
was selected to be near the afterglow position, but located so that it did not contain any obviously detected stars. 
Source counts were determined by subtracting the expected average background counts inside the source extraction 
region from the total counts inside that region. The v band count rates are shown as a function of time since burst in 
the top panel of Fig.[6j The best fitting power law was determined using the technique described at the end of § ??. 



The uncertainty on the power law slope was estimated by calculating the rms of the best slope values in 100 simulated 
data sets. In each of the simulations the counts in the source and background regions were randomly varied according 
to a Poisson distribution, with means equal to the measured values. The best power law slope and its uncertainty are 
a = —0.62 ± 0.19; the best fitting power law is shown i n Fig. [6] The best slop e value is consistent with the value found 
from the UVOT whit e band data of a = —0.62 ±0 . 32 (iBlustin fc Page 2006), and from the ground-based R band data 
of a = -0.35 ± 0.36 (jTanvir. Rol. fe rlewetdl2006l : lGarnavichll2006h . It is clear that although none of the individual 
measurements give a significant detection, the power law decay function fit to all of the data is quite reasonable. In 
this case, we have used the fitted value of a because it is consistent with the decay slopes in other bands for which 
there are afterglow detections, and by design it will give the maximum weighted S/N. In the absence of a reliable fit 
in the v band, it would have been appropriate to adopt a value of a given by the results in the other bands. 

The cumulative S/N at each exposure was calculated from equation ([4]) for the unweighted data, and from equation 
(|23|) for the weighted data. The bottom panel of Fig.[6]shows the cumulative S/N for both the weighted and unweighted 
cases. The S/N rises with coaddition in both cases up to the first few images. In the unweighted case the S/N drops 
from a high of about 2.8 to about 2.0 after coadding all of the data. In the optimal weighting case, the S/N rises to 
about 3.0, and remains quite steady ther eafter. Both the maxim um unweighted S/N and the final weighted S/N are 
greater than the value (2.4) reported bv lBlustin fc Pa ge (2006). Using optimal coaddition, the v band data by itself 
would have been sufficient to detect the optical afterglow. The individual images were coadded using no weighting, 
and using the optimal weights calculated as described here. Images centered on the afterglow position are shown for 
both cases in Fig|7l The images have been scaled so that the field stars appear to have the same brightnesses in both 
images; this allows a direct comparison of the relative coadded brightness of the afterglow in both images. The effect 
of optimal coaddition is to increase the S/N of the afterglow relative to unweighted coaddition. The process also 
degrades the S/N of other objects in the image, since for them the weights are not optimal. The afterglow is difficult 
to see in the unweighted coadded image, but easily seen in the weighted coadded image. 

This case illustrates the potential of optimal coaddition to improve the dete ction rates of GRB afterglows. We plan 
to apply this method to the entire UVOT GRB database (Ro ming et al.|[2008l ). to determine if new detections of other 
afterglows can be achieved. 

5.3. Improved Detection Limits for Optical Non-Detections 
Op tical afterglows have not been detected for many GRBs to date, including GRBs discovered by Swift (|Roming et al.l 



2006). It is important for detection limits to be as deep as possible in order to help understand the "dark burst" 
phenomenon. Also, when attempting to identify very high redshift burst candidates, fainter detection limits in bluer 
bands place tighter constraints on the range of possible redshifts. Here we examine the case of GRB 060923A, which 
is a potentially very high redshift burst, with detections in the K band, but no detections at optical bands, including 
the UVOT v band. 

GRB 060923 A was discov ered by the Swift BAT , and the UVOT began observations of the burst location about 
85 seconds after the trigger (jStamatikos et al1l2006f) . No afterglow was detected in any of the UVO T images, either 
singly or coadded, and no afterglow was reported from other observations in the V, R, I, and J bands (|Li et al.ll2006bl : 



iMelandri fc Gombod [20061: 1 Williams fc Milne) 120061 iFox. Rail, fc Of ck 2006). An aft erglow was found and confirmed 
only in the K band (jTanvir. Rol. fc Hewettl 120061 : IFox. Rau. fc Ofekl 120061 : IFox! 12006). suggesting that GRB 060923 A 



could be at a very high (z > 10) r edshift or heavily reddened. Subsequent very deep R band observations revealed a 
possible host galaxy with R ~ 25.5 (jTanvir. Rol. fc rlew ctt 2006), which would favor the reddening hypothesis. Tighter 
constraints on the bluer band magnitudes would provide correspondingly tighter constraints on either the redshift, or 
extinction of the GRB. 

As discussed in § 13.31 optimal coaddition can be used to find the maximum count rate (and therefore deepest 
magnitude) that an afterglow could have had in the first image (where the count rate would have presumably been the 
greatest) without being detected at a given significance level. The detection limit calculations require the exposure 
start and stop times, an estimate of the afterglow decay rate, and a measurement of the background counts at the 
GRB position in each image. UVOT images of the burst position were taken in the v band from about 85 s to about 
11 ks after the trigger. Background counts in each of the images were estimated for a three arcsecond radius aperture 
at the reported location of the K band afterglow. The temporal decay rate is unknown, but reasonable estimates can 
be made from several arguments. First, the K magnitude of the afterglow dropped by at least 0.9 magnitudes from 
about 6.7ks to 86 ks after the burst, which means the K band decay slope was less than an = —0.33. Second, the 
XRT li ghtcurve faded at two distinct rates: cvxrti = —2.7 from about 80 to 400 s, and oixrti = —1-23 after about 
4000 s (jConciatord I2006D . These values bracket typical optical afterglow decay slopes; calculations were made using 
all three values. 

Detection limits for a S/N of three in the first UVOT v band image were made for the image by itself, the unweighted 
coadded image, and the optimally weighted coadded image, as described in § 13.31 The count rate li mits were converte d 
to standard UVOT v band magnitudes using the zero points and aperture corrections given bv iPoole et all (|2008h . 
The results for each of the three decay slopes are given in TablejTJ The single image limit is v — 20.22. For unweighted 
coaddition, there is little (v = 20.26) or no improvement with a shallow decay slope (a = —0.33), and the limit 
actually becomes worse with steep decay slopes. The reason for this is that noise in the later images degrades the S/N 
of unweighted coadded images, so that a brighter source in the first image is required to produce a significant detection 
in the final coadded image. For a shallow slope, the optimally weighted v band limit (v — 20.42) in the first frame 
is significantly deeper than the single frame and unweighted coaddition limits. For steep decay slopes, the optimally 
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weighted limit is virtually the same as the single frame limit (and much better than the unweighted coaddition limit). 
There is little improvement over the single frame limit in this case, because the steep lightcurve decay makes the later 
image weights very small, so that they contribute very little to the final S/N. 

The early UVOT v band detection limits are the deepest reported for GRB 060923A. Any estimate of either a 
photometric redshift, or an extinction value for this burst, must be consistent with the detection limits. The reported 
K band detection measurement was made more than 6200 s after the first v band measurement. Given the range 
of decay slopes considered, the v band magnitude limit at the time of the K band measurement would be at least 
1.3 magnitudes fainter than in the first image. This would place even tighter constraints on either the redshift or 
extinction. The same analysis can be applied to time series images of any non-detected GRB afterglow to improve the 
detection limits. 

6. DISCUSSION AND SUMMARY 

Image addition for the detection of variable objects, or the improvement of S/N , can benefit significantly from 
proper weighting. Data from rapidly varying sources, such as GRB afterglows, stands to benefit the most from 
optimal coaddition. In fact, it can be counterproductive not to optimally weight the data when coadding. We have 
shown how optimal coaddition can ensure the best S/N for rapidly varying GRB afterglow data, even to the point of 
recovering previously undetected sources, and can place tighter detection limits on image data sets. 

We have assumed a s imple power law model for the optical lightcurve; this will clearly be inadequate in many 
cases (|Oates et al.l [2008). A spectral extrapolation from the X-ray to the optical lightcurve may provide a reasonable 
estimate at early times before the for ward shock peak . The spectral index (3 should vary between (3 — (p — l)/2 and 
(3 = p/2 where p is the electron index (|Sari et~allfl99l . Given a range of possible p values, one could predict a range 
of possible optical lightcurves. However, at early times in many X-ray lightcurves, there is a fast decay phase with 
superimposed flaring activity, which rarely has a c ounterpart at opti cal wavelengths. It is probably better in almost all 
cases to assume an average optical lightcurve fe.g. lOates et al.ll2008l ) rather than extrapolate from an X-ray lightcurve. 
In any case, it is straightforward to adapt the coaddition technique to arbitrary lightcurves. 

The technique can also be applied to observations of any variable source (e.g. supernovae, active galactic nuclei, 
cataclysmic variable stars), as long as a reasonable estimate of the lightcurve is available. Optimal coaddition is 
most beneficial when count rates are low, background is high, and the rate of source variability is high. All of these 
factors often apply to optical observations of GRB afterglows. Optical observations have been used in the examples 
presented here, but the technique is equally applicable to UV, infra-red, and X-ray imaging data that has Poisson 
noise characteristics. The benefit to X-ray and possibly UV data may not be as great as for other frequency ranges, 
since the background count rates are often quite low. However, there would likely be some improvement in S/N or 
detection limits, and no loss in either, when the source count rates are comparable to the background rates. 

There are other ways of improving S/N that could probably be used in combination with the optimal coaddition 
technique described here. For example, we have used simple circular apertures to extract source counts in individual 
images; if the image point spread function can be determined with precision, there is the potential for optimally 
extracting the cou nts in each im age, or the final coadded image. A related optimal pixel weighting technique has 
been described by iNavloil |l998). In addition, if the characteristics (aside from source brightness) of images differ 
significantly from frame to frame (which is usually not the case f or UVOT data), there are techniques to optimally 
coadd the images to achieve improved S/N or image resolution (jFischer fc Koc hanski 1994). As another example, 
the UVOT can observe in event mode, by tagging the arrival time of individual detected photons. This allows the 
possibility of optimally weighting individual events, rather than integrated images, which could improve the S/N of 
afterglow data, particularly at early post-burst times. We also caution that the equations derived here assume count 
rates high enough that Gaussian noise statistics can be assumed, but at extremely low count rates that assumption 
may not apply, and other approaches should be considered. 

We will apply optimal image coaddition to all of the GRB afterglow data obtained with the Swift UVOT. The 
primary goal will be the recovery of previously undetected afterglows, but at a minimum, the detection limits will be 
improved, placing stronger constraints on the existence of so-called dark GRBs. 

Funding for the Swift program at Penn State is provided by NASA under the contract NAS5-00136. We acknowledge 
the use of public data from the Swift data archive. 
Facilities: Swift (UVOT) 
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TABLE 1 

GRB 060923A UVOT v band magnitude limits for the first image. 
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Fig. 1. — Cumulative S/N as a function of time since the burst, for coaddcd simulated GRB afterglow data, with a shallow decay slope 
(left) and a steep decay slope (right). The temporal decay slopes a, and the initial count rate Hi, of the lightcurve models are given, along 
with the background count rate R a k y - Each point corresponds to an observation in a relatively short time interval. Unweighted coadded 
data is shown by red triangles; optimally weighted coadded data is shown by blue squares. The error bars span the 68.3% confidence limits 
in the results of 1000 simulations. 
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Fig. 2. — Cumulative S/N as a function of time since the burst, for coadded simulated GRB afterglow data, with a 
initial count rate (left) and high initial count rate (right). The temporal decay slope a, and the initial count rates of 
models are given, along with the background count rate R a k y - Each point corresponds to an observation in a relatively short 
Unweighted coadded data is shown by red triangles; optimally weighted coadded data is shown by blue squares. The error 
68.3% confidence limits in the results of 1000 simulations. 
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Fig. 3. — Cumulative S/N as a function of time since the burst, for coadded simulated GRB afterglow data, with a high background 
count rate (left) and a low background count rate (right). The temporal decay slope a, and the initial count rate Ri, of the lightcurve 
models are given, along with the background count rates R s ky Each point corresponds to an observation in a relatively short time interval. 
Unweighted coadded data is shown by red triangles; optimally weighted coadded data is shown by blue squares. The error bars span the 
68.3% confidence limits in the results of 1000 simulations. 
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Fig. 4. — Final S/N in the optimally coaddcd data as a function of the model temporal decay index a mo£ j e i, for simulated GRB afterglow 
data. The true decay index is —1.2, and the simulations used an initial count rate of 3ctss — 1 , and a background count rate of 2.5ctss~ 1 . 
The points give the median values from 1000 simulations, and the error bars show the uncertainty on the median value (not the dispersion 
in the simulations). The dashed line shows the maximum S/N from unweighted coadded data averaged over the simulations. 
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Fig. 5. — Lightcurve (top), and cumulative S/N (bottom) for UVOT v band image data of GRB 050525A, as a function of time since 
the burst. Error bars on the lightcurve represent count rate uncertainties for measurements using individual images. In the S/N plot, red 
x's show the results of unweighted coaddition; blue squares show the results of optimal weighting coaddition. 
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Fig. 6. — Light curve (top), and cumulative S/N (bottom) for UVOT v band image data of GRB 060604, as a function of time since the 
burst. Error bars on the lightcurve represent count rate uncertainties for measurements using individual images. In the S/N plot, red x's 
show the results of unweighted coaddition; blue squares show the results of optimal weighting coaddition. 





Fig. 7. — Coadded UVOT v band images of GRB 060604. Left — result from unweighted coaddition. Right — result from optimally 
weighted coaddition. Each image is 4 X 4 arcmin on a side. The images have been scaled so that the brightnesses of the field stars appear 
nearly the same in both images. The reported location of the afterglow, detected in other bands, is at the center in each case. The circles 
show the 3arcsec radius apertures used to extract the source counts. In the unweighted case, the S/N of the counts in the source region is 
about 2.0, while in the weighted case it is about 3.0. 



